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Abstract 

Starting with a finite A;- mesh version of a well-known equation by Blount, we show 
how various definitions proposed for the polarization of long chains are related. Ex- 
pressions used for infinite periodic chains in the 'modern theory of polarization' are 
thereby obtained along with a new single particle formulation. Separate intracellular 
and intercellular contributions to the polarization are identified and, in application 
to infinite chains, the traditional sawtooth definition is found to be missing the 
latter. For a finite open chain the dipole moment depends upon how the chain is 
terminated, but the intracellular and intercellular polarization do not. All of these 
results are illustrated through calculations with a simple Hiickel-like model. 
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The purpose of this Letter is to answer some questions that arise in connec- 
tion with the theoretical treatment of macroscopic polarization in quasi-one 
dimensional chains. In order to specify the issues let us consider a macroscopic, 
open-ended, polymeric chain consisting of identical unit cells. For sake of ar- 
gument this chain is assumed to be polarized due to an asymmetric unit cell 
and/or an external electric field. There are two contributions to the polariza- 
tion P, i.e. to the dipole moment per unit cell. One is due to the asymmetric 
charge distribution within a unit cell in the central region of the chain and 
the other is due to the charge of opposite sign that accumulates at the chain 
ends. Hiickel-type calculations (see below) show that the contribution due to 
the finite chain ends does not vanish even in the infinite chain limit. This is 
simply due to the fact that, for a charge of fixed magnitude at either end, the 
dipole moment is directly proportional to the distance between the charges. 
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Next, imagine that the chain ends are connected to form a ring. In that event 
there are no ends and all the unit cells are identical. What is the relation- 
ship between the unit cell charge distribution of the closed chain and the unit 
cell charge distribution at the center of the open chain? In fact, they are the 
same as our Hiickel-type calculations confirm. Then, what has happened to 
the contribution to the polarization associated with the charge build-up at the 
ends of the open chain? As it turns out this contribution is associated with 
a charge flow term that arises from Blount's theoretical expression [1] for the 
polarization when periodic boundary conditions are applied. This raises the 
question of whether or not such a term can be accounted for by the conven- 
tional sawtooth approach [2,3]. The latter is based on using a finite mesh in 
/c-space, along with periodic boundary conditions [4,5], but it does not corre- 
spond to a finite-mesh analogue of Blount's formula which will be presented 
here. From a general formulation of this analogue several approximations will 
be developed including the fundamental equation(s) of the so-called modern 
theory of polarization [6,7,8,9,10] 

Using a Hiickel model, discussed below, we demonstrate quantitatively that 
the sawtooth approach omits the current term and that Blount's formula gives 
an accurate approximation for the polarization when applied with a finite set 
of k points. Finally, it will be seen that the polarization of long finite chains 
with arbitrary terminal substituents does not depend on the nature of these 
substituents even though the same cannot be said of the dipole moment D 
itself. 

For a finite chain of alternating A and B atoms the Hiickel Hamiltonian may 
be written in terms of orthonormal atom-centered basis functions {Xp} as 

2K 2K-1 

H= ^ e p£\pp~ *p,p+i(Vh^p + c£c p+ i) (1) 

p=-2K+l p=-2K 

where 4K is the number of atoms, c P and c p are the creation and annihila- 
tion operators for the function Xp, an d e p and — t PtP+1 are on-site energies and 
hopping integrals, respectively. The matrix elements of the position operator 
are given by {Xn\z\Xm) = ^n,m z n with z p being the position of the pth atom. 
For convenience, the atoms are taken to be equally spaced, z p = |[p — |], 
with p = —IK + 1, . . . , 2K being odd (even) for the A (B) atoms, and | 
being the nearest-neighbor distance. Assuming one electron per atom it is 
straightforward to evaluate the atomic charges and the electronic polariza- 
tion, P (dipole moment per A-B unit) as a function of the number of A-B 
units, 2K. As an example (using arbitrary units for length and energy and 
setting the electronic charge equal to +1) for a = 2.0, alternating on-site en- 
ergies e p = ±0.5 = ±e , and hopping integrals t + = —2.2, i_ = —1.8, we 
find that the polarization is converged to a value P = 0.58125 for K > 13 
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while the charges on the central atoms are Q(A) = 0.74413, Q(B) = 1.25587. 
The latter result in a contribution P c = P - f [Q(B) - Q(A)] = 0.32538 to 
the polarization due to accumulation of charge at the chain ends. We may 
think of P c as an intercellular charge flow term. The existence of a substantial 
intercellular charge flow term is remarkably robust to variations in the model. 
Thus, including next-nearest-neighbor interactions, modifying the matrix ele- 
ments at the chain ends or altering (Xn\z\Xm) in realistic ways often changes 
the total dipole moment, but not the polarization, of sufficiently long chains. 
The effect of varying terminal (on-site and/or hopping) matrix elements is 
illustrated in Fig. 1. The figure shows that altering the chain ends changes 
the charges on the ends (upper part) and the dipole moment (lower part) but 
neither the polarization (slope of dipole-moment curve) nor the charges in the 
central region are affected. The right hand panel of the lower part shows that 
the polarization can vary for small chains (see, particularly, curve f). 

The fact that one may calculate the polarization by studying only the central 
cells has been shown previously by Vanderbilt and King-Smith [9]. However, 
this does not imply that different terminating groups, which lead to different 
charge accumulation at the chain ends, give the same polarization because 
there is a contribution due to intercellular charge flow that could change. 
Thus, this is a generalization of the Vanderbilt and King-Smith result, which 
is consistent with the known near-sightedness [11] of the single-particle density 
matrix. It has obvious implications for the design of donor-acceptor, or push- 
pull, systems and is valid provided the chain is sufficiently long. 

Even for an unsubstituted chain where the ends are connected so that no 
charge can accumulate there is an important contribution to the polarization 
that arises from the intercellular charge flow. We now turn to that case and 
consider a ring of IK identical AB unit cells. Application of periodic boundary 
conditions leads to the general expression for the eigenfunctions 



with n being a band index, Xpm the pth basis function of the rath unit cell, and 
Nb the number of basis functions per cell. In our Hiickel model, n = 1 (n — 2) 
for the occupied (empty) band, and p = 1 {p — 2) indicates the function on 
the A (B) atom. For any given set of parameters and sufficiently long chains 
the electronic charges, Q(A) and Q(B), turn out to be identical to those at the 
center of the open-ended chain of the same length. Thus, for either chain, the 
same intracellular polarization |[Q(B) — Q(A)} is obtained. This means that 
the P c contribution must be accounted for in some other way. For the ring we 
can readily identify that contribution by considering Blount's expression for 



ikz 




(2) 
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the polarization in the limit K — > oo: 

^f?/^!"** (3) 

where the n summation is over the (doubly) occupied bands. Using Eq. (2) it 
is easy to show (cf. [12]) that 



l r ( d \ 

P=~J2 E / e ' fema E c ?n ({x g o\z-ma\x pm )+t(x q o\Xp m )^r )4n dk - 

71 n m =-N J PQ V a/C/ 

(4) 

With a finite 2-fT-point-mesh approximation for the integral one can verify that 
the first term on the rhs of Eq. (4) yields the intracellular polarization. This 
leaves the second term as the periodic cyclic chain analogue of the intercellular 
charge flow contribution described above in connection with the open chain. 

If one is interested in long open-ended chains it is usually advantageous com- 
putationally to assume that the chain is infinite and periodic. A number of 
different proposals have been advanced for calculating the polarization of in- 
finite periodic chains using finite k mesh methods. In order to compare these 
approaches we follow the treatment of Blount [1], based on the relation 

ztfjf) = ie ikz ^-e- ikz ^ - * J-^(r-), (5) 



to obtain the effect of the coordinate z acting on a single electron whose 
orbital, ip(r), is expanded in terms of Bloch waves 

m = EE^U* = EE^n(^/« fc (6) 



In fact, Blount [1] obtained Eq. (4) by using Eq. (5) on ip of Eq. (6). Here 
instead of a continuous k we will use a finite /c-mesh, which corresponds to 
assuming that the system possesses the periodicity of the Born von Karman 
(BvK) zone containing 2K unit cells. Consequently, the analytical derivatives 
in Blount's formulation will be replaced by numerical derivatives. In lowest 
order the numerical derivatives corresponding to the terms in Eq. (5) are: 
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A^ = ±(A'_ + A' + ) 
A» = i(A»+A»), 



(7) 



with Ak = By construction these expressions have the BvK periodicity 
and it follows that the lowest order finite- k- mesh analogues of (ip\z\ip) are 



^IHAV + .A';)|^) = (^| — (e 



^|(-<Ai + tAg)|^> 



^-1)1^) : 

sin(AA;z) 1 



(8) 



where S ± = {^\e ±lAkz \^). If the spatial extent of ip is much smaller than 
(this can, e.g., be obtained by increasing the number of k points in an 
actual calculation), and assuming that J BvK \zip{f)\ 2 df exists (e.g. when ip is 
a well-localized Wannier function), then we may make the approximation 
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(9) 



or 



ty\ z \f)~-± l n S + ~ -Um\nS + 

ZA/C ZA/C 

i _ i 

(i/j\z\ip) ~ — — ln S ~— — ImlnS* 

ZA/C ZA/C 

(V>Nv> ^ ^ Arcsin t^( 5+ - s-)]. (io) 

Here, the second equalities in the first two expressions have been obtained 
by removing the imaginary parts and, accordingly, requiring that (ip\z\ip) is 
real. This result comes about automatically in the expressions based on the 
A operators. 

The treatment for N electrons is similar. In that case the one-electron Bloch 
waves are replaced by Slater determinants ^ k = A[ip k ^ (ri)^ 2 (r 2 ) • • • ip k ^ (f/v)], 
where A is the antisymmetrizer, and Eq. (5) becomes 
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Then an arbitrary iV-electron function can be written as the linear combi- 
nation ^/(rl, r 2 , . . . , r N ) = X^X^ ^(^1^2, ■ ■ ■ , r N )f£ with the single-particle 
situation being a special case. The generalization of the quantities in Eq. (7) 
becomes 



A'_*(ri,f2,...,rW) 



A'LWi,r2,...,r N ) = 
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with analogous expressions for A' + and A" , A' and Aq. Hence, the general- 
ization of Eq. (8) is 



A_ = -iA'_ + iA"_ 
A+ = -iA' + + iA" + 

A = -iAg + iAg' = 



I 

~~ Ak 
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_ e iAfc(zi+^ 2 H hzjv) 
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Ak(z 1 + z 2 -\ h zjv) 



(13) 



We will restrict ourselves to the case where there is a finite gap between 
occupied and unoccupied bands, assume no spin polarization (N is even), and 
use a single determinant wavefunction (Hartree-Fock or Kohn-Sham theory). 
Then 



(tf|A_|tf) = 


Ak 


1- (det£+) 2 
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Ak 


(detif) 2 - 1 



(*|Ao|*> 



HAk 



(det£ + ) 2 -(det£-) 2 



(14) 
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where S_ is the N/2 x N/2 matrix containing the single-particle matrix el- 
ements St k \ un = (il)i\e ±lAkz \il} l j) . Assuming localized orbitals we may apply 
the analogue of Eq. (9), i.e. 

/^yi 1 ±iAk(zi+Z2+-+z N )\q f \ ^ 1 ±iAk(S>\zi+Z2+-+z N \y) r\ 

Ak 1 } ~Ak ' 1 ' 



and either of the first two equations in Eq. (14), in combination with Eq. (13), 
to arrive at the expression 

P R = --Im In detS + = -Im In det£T, (16) 

7T — 7T 



for the polarization. Note that {S_ + y = S.~ ■ An essentially identical formula in 
terms of Bloch orbitals has been given by Resta [4]. We observe that Eq. (16) is 
based on a not too accurate finite-difference approximation to the derivative. 
A more accurate approximation is 



— Arcsin 
2tt 



l((detr) 2 "(detr) 2 ) 



(17) 



Despite this Pr turns out to be more useful computationally. To see why we 
write 

det^ = s±it, (18) 
whereby 



P R = Arctan(-) 

7T S 

P = — Arcsin(2st). (19) 
2tt 

As K — > oo, s 2 + t 2 — > 1, while \s\,\t\ < 1 are increasing functions of K. 
Accordingly, as our numerical results below confirm, Pr converges faster than 
P as a function of K. On the other hand, P can be valuable analytically; 
indeed, it motivated our choice for the operator defined in Eq. (23) below. 

The value of det Sj^ will not be altered by an arbitrary unitary transformation 
of the single determinant orbitals. So, instead of localized orbitals we may use 
the occupied Bloch waves from which these orbitals are obtained. Then, the 
matrix elements of Sj^ are non-zero only for pairs of Bloch waves whose k 
values differ by Ak (modulus ^f). As a result ^ ± can be written as consisting 
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of 2K x 2K square blocks, each of dimension B — (the number of doubly- 
occupied bands) with non-zero elements only in the set of blocks lying one 
stripe above and one stripe below the main diagonal. 

Given that the Bloch functions are differentiable with respect to k as discussed 
by Blount we obtain for small Ak 

2K B o 

ln(det^) 2 ~ T 2AkJ2 E (20) 

fc = l 71=1 



Inserting this into the rhs of Eq. (16) yields another formula for the polariza- 
tion 

• 2K B o 

^ksv = ^EE(^I^), (21) 

fc=l n=l 



which is the ID discretized Berry phase expression [9] used in the modern 
theory of polarization. 

For the treatment of core orbitals (or those of non-interacting periodically 
repeated molecules) we suppose that the orbitals are strongly localized so 
that {^> Pl (r — R ni )\e ±lAkz \ip P2 (r — R n2 )) vanishes unless the units n\ and n 2 , 
where the functions are centered, are identical. In that case we may write for 
orbitals of the same unit 



(VV \e ±iAkz \*P P2 ) ~ S puP2 e ±AkZp ^ ± iAk(4> Pl \z - z puP2 \i> P2 ) 

+ W™Lty pi \( z - z PuP2 ) 2 \^ P2 ) + ■■■ (22) 

with z PljP2 being the 'center' of the pith and p2th orbital. Therefore, the 'tradi- 
tional' contribution to the polarization from these orbitals, i.e., Sp^pMV'p)) 
is obtained only in the case where all terms but the first one on the rhs of Eq. 
(22) are negligible (e.g. in the limit AA; — > 0). 

So far we have presented an internally consistent approach for how to calculate 
the polarization in an infinite, periodic chain when basing the discussion on a 
generalization of Blount's work to the case of a finite BvK zone. We have ar- 
rived at an expression involving the expectation values for iV-body operators, 
i.e., the Sj^ matrices. This has been taken as a proof that the polarization is a 
many-body phenomenon [6,7]. However, the polarization can also be written 
in terms of the single-particle operator 
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It is straightforward to show that the expectation value of this operator gives 
P in the limit AA; -> 0. 

In order to explore our ideas further, the Hiickel-like model described above 
was used to evaluate the various polarization expressions we have presented. 
For our purposes it is necessary to have matrix elements of z and e ±tAkz that 
are defined consistently. Hence, we calculated the matrix elements in both 
cases analytically assuming piecewise constant basis functions of adjustable 
width, w. For simplicity we also assumed that w < |, whereby the results be- 
come independent of w. Other more realistic functions are possible, of course, 
but the above choice is sufficient to make the desired comparisons. 

In Table 1 we show some typical results for the various choices of P obtained 
using BvK periodic boundary conditions. The finite chain value determined 
from the increment AD = \[D(2K + 2) - D(2K)], where D is the dipole 
moment, is also presented for comparison. In order to interpret polarization 
values the reader should recall that P is determined only up to an arbitrary 
multiple of the unit cell length [cf. Eqs. (16) and (17)], which in this case is 
2.0. Bearing this in mind, the table shows that AD agrees very well with the 
polarization of the infinite system given by Pr. Indeed, the finite chain result 
converges more rapidly to the infinite K limit. The sawtooth approximation 
(denoted P st in the table) is calculated using periodic boundary conditions 
with z replaced by a piecewise linear function having the BvK periodicity. 
Note that P st gives the correct value only when the system consists of purely 
non-interacting units (last case in table). Since sin(a) = sin(7r — a) it is not 
possible to discriminate between P and | — P (cf. the first case in the table). 
If that is taken into account, we see that -Pksv> Pr, an d Pq all give similar 
results, although the latter converges much slower, and the former much faster, 
than the others. 



In conclusion we have provided a unified picture of electronic polarization in 
extended quasilinear chains based primarily on the finite /c-mesh analogue of 
Blount's treatment for infinite periodic systems. Separate intracellular and 
intercellular contributions are identified and compared between closed and 
open chains. It is shown that neither component is affected by substitution at 
the end of an open chain, as occurs in a push-pull compound. On the other 
hand, the traditional sawtooth formulation for infinite closed chains fails to 
account for the intercellular charge flow term. Several different expressions for 
the electronic polarization are systematically generated from the same starting 
point, including those related to the so-called modern theory of polarization. 
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From the same perspective we obtain an alternative single particle operator, 
which yields the polarization as its expectation value. Hiickel-type calculations 
are carried out to illustrate all of these points and to assess the convergence 
properties of the various polarization formulas as the /c-mesh spacing decreases 
to zero. 

This work was supported by the German Research Council (DFG) through 
project Sp 439/11. One of the authors (MS) is grateful to Fonds der Chemis- 
chen Industrie for generous support. 
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Table 1 

Results of model calculations with the Hiickel model. The lattice constant equals 
a = 2. All other parameter values are given in the table. 





t+ 


t- 


K 




Pr 


-Prsv 


^0 




AD 


0.5 


2.2 


1.8 


20 


0.25587 


-1.41859 


-1.41745 


0.29410 


0. 


.58125 








200 


0.25587 


-1.41875 


-1.41875 


0.39798 












2000 


0.25587 


-1.41875 


-1.41875 


0.41643 




» 








20000 


0.25587 


-1.41875 


-1.41875 


0.41852 






0.5 


2.5 


1.5 


20 


0.21337 


-1.68291 


-1.68305 


0.26857 


0. 


.31695 








200 


0.21337 


-1.68305 


-1.68305 


0.31134 




» 








2000 


0.21337 


-1.68305 


-1.68305 


0.31638 












20000 


0.21337 


-1.68305 


-1.68305 


0.31690 
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1.5 


1.5 


20 


0.33562 


-1.00000 


-3.00000 


0.00000 


1, 


.00000 








200 


0.33562 


-1.00000 


-3.00000 


0.00000 
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2000 


0.33562 


-1.00000 


-3.00000 


0.00000 




J? 








20000 


0.33562 


-1.00000 


-3.00000 


0.00000 
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1.5 


20 


0.00000 


-2.00000 


-2.00000 


0.00000 


0. 
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200 


0.00000 


-2.00000 


-2.00000 


0.00000 












2000 


0.00000 


-2.00000 


-2.00000 


0.00000 












20000 


0.00000 


-2.00000 


-2.00000 


0.00000 






0.5 


2.0 


0.0 


20 


0.24254 


-1.75735 


-1.75746 


0.22586 


0. 


.24254 








200 


0.24254 


-1.75746 


-1.75746 


0.24078 
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2000 


0.24254 


-1.75746 


-1.75746 


0.24236 




?? 








20000 


0.24254 


-1.75746 


-1.75746 


0.24252 
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Fig. 1. Upper part: Distribution of the atomic charges as a function of atom index 
for a finite chain with 80 atoms for different cases of the matrix elements for the 
terminating atoms, i.e., the on-site energies for the first (last) atom have been mod- 
ified as eo — > eo + Aei (— eo — ► — eo + Ae r ), and the hopping integrals between the 
first (last) two atoms according to — 1+ — ► — 1 + + Ati {—t + — ► — 1+ + At r ). Lower 
part: The dipole moment for the same cases but as a function of chain length. 
Here, 'sawtooth' corresponds to the dipole moment for a ring system when using 
the sawtooth approximation, and the right panel shows a blow-up of the low-K 
part. The curves marked a, b, c, d, e, and / (in the lower panels these labels are 
listed in the same order as the curves appear) correspond to the following modi- 
fications: (Aei,Ae r ,Ati,At r ) = (0,0,0,0), (1.0,0,0,0), (0,1.0,0,0), (0,0,-1.0,0), 
(0,0,0,-1.0), and (1.0,-1.0,-1.0,1.0), respectively. 
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